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DEVICE FOR PRODUCING DYNAMIC DISTORTION PATTERNS 
AT INLETS OF AIR-BREATHING ENGINES 
by Robert J. Baumbick 
Lewis Research Center 

SUMMARY 

The development of the air distortion system is discussed. Also included is a dis- 
cussion of the servosystem designed to drive the air valve. An analysis of the pneumatic 
system dynamics was made and correlated to experimental results. The objective in 
the design of the air distortion system was to produce pressure distortions at the engine 
face. This is brought about through a momentum exchange between primary air flow to 
the engine and secondary flow produced by the distortion valve. Problem areas that 
affect jet flow are cited, and corrections are suggested for improving the response. 

Both analytical and experimental results showed how supply line dynamics degrade the 
response of the system. From these results, future systems will be modified to reduce 
these supply line dynamics. Experimental and analytical results compared quite well 
for both the servosystem and the pneumatic system. 


INTRODUCTION 

The trend in engine development is toward increased performance without substantial 
increase in engine weight. Weight savings are realized by increasing the work output 
per stage of the compressors. The higher stage loading of these compressors results in 
a greater sensitivity to inlet flow distortions. Information on the effects of inlet distor- 
tion on engine performance is presented in references 1 to 3. 

Because of the higher stage loading, the compressor is forced to operate closer to 
its stall region. Research engineers are concerned with defining the engine stall region 
when flow distortions are introduced to the inlet of the engine. In the past, uniform and 
nonuniform pressure disturbances have been produced by placing screens in the engine 
inlet duct. Cyclic distortions were produced with devices such as the rotating bypass 
valve (ref. 4) and the multispoked rotating valve (ref. 5). These devices, however, were 



capable of producing only sinusoidal pressure disturbances. Since the exact nature of 
inlet disturbances that affect engine performance are not well defined, a more versatile 
tool was needed for engine, distortion studies. This need led to the development of the 
air distortion system discussed in this report. 

Six valves were installed at 60° intervals around the engine inlet duct. Figure 1(a) 
shows the air valves installed in the engine test chamber. Each valve discharges into 
nine tubes. These tubes are terminated in orifices. The discharge jets (54 total) form 
the array shown in figure 1(b) as it appears looking toward the engine inlet. Indicated 
in this figure is a sector of nine discharge jets. The system schematic, together with 
the location of the discharge jet plane and the pressure measuring station in the inlet 
duct, is shown in figure 1(c). 

Pressure distortions are produced at the engine inlet through a momentum exchange 
between primary engine flow and secondary flow produced by the distortion system. Be- 
cause each valve can be driven independently, a wide variety of pressure distortions can 
be produced at the engine inlet. These distortion patterns can be steady state, sinusoi- 
dal, or transient in nature. Wenzel (ref. 6) presents results of transient pressure dis- 
tortions produced at the engine inlet by driving the air valves with discrete pulses. 

These pulses varied in both amplitude and duration. Circumferential distortion patterns 
from 60° to 360° in 60° sectors were produced at the engine inlet. 

Meyer (ref. 7) presents results of steady -state and sinusoidal pressure distortions 
produced at the engine inlet. The pressure amplitudes produced at the engine inlet were 
attenuated beyond 20 hertz. The poor distortion patterns produced at the engine inlet at 
higher frequencies were the result of inlet duct dynamics and distortion system dynam- 
ics. The air distortion system performance is the subject of this report. Evaluation of 
the air valve and discharge jets was accomplished in bench tests using a prototype of the 
distortion system used in the engine test facilities. 


APPARATUS AND PROCEDURE 

The prototype test rig used for evaluation of the distortion system is shown in fig- 
ure 2. It consists of an accumulator, supply line, air valve, and nine discharge tubes 
that are terminated in orifices. Pressure measurements are made at the inlet and dis- 
charge ports of the valve and at a point upstream of the discharge tube terminating ori- 
fices. 

The air valve, shown in figure 3, consists of two concentric cylinders, each having 

36 rectangular slots cut out on their circumference. Valve area is a linear function of 

the relative position between the two cylinders. A maximum valve stroke of 0. 25 inch 

o 

(0. 635 cm) yields a flow area of 7. 96 square inches (51. 3 cm ) and a flow of approxi- 
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mately 7 pounds per second (3. 175 kg/sec). The valve is driven by a fast -response (150 
to 200 Hz) two -stage electrohydraulic servovalve assembly. The servosystem was de- 
signed at the NASA Lewis Research Center and incorporates a number of unique features 
that extend the bandwidth of operation of the basic servovalve. More information on the 
design of the servosystem is presented in reference 8. The air valve position is sensed 
by an LVDT (linear variable differential transformer) and is used as the feedback signal 
for closed-loop control of the air valve. 

Unbonded strain-gage pressure transducers were used for pressure measurements 
at the valve inlet and discharge. These were the differential pressure type and were ref- 
erenced to the atmosphere. The transducers have a frequency response of better than 
5 kilocycles with accuracy of better than 1 percent. A quartz -type pressure transducer 
was used for pressure measurements at the discharge orifice. These transducers with 
connecting cables have a frequency response of better than 5 kilocycles. The transducer 
used to measure the valve position has no dynamics in the frequency range of interest.. 

Steady-state data at a given operating point were taken on strip-chart recorders. 
Frequency responses were run by using a sweep-frequency method. A sweep rate of 
1 decade per minute was used. The time-varying portions of the signals were recorded 
on tape. The data were reduced by the method described in reference 9, which also 
discusses the errors incurred in using the sweep -frequency technique. 


Distortion Valve Position Servo 

A block diagram of the position servo used to drive the air valve is shown in figure 4. 
The major components are the servoamplifier, the servovalve, and the load, which con- 
sists of the piston-in-cylinder actuator and the air valve. The first stage of the servo- 
valve is an electrical torque motor that converts an input current into a torque. The 
torque actuates a flapper valve that, in turn, drives a spool valve. The spool valve 
directs high-pressure hydraulic fluid to the load actuator. The dynamics of the torque 
motor coils, due to the inductance and resistance of the coil, are included in the servo- 
amplifier section. The second block contains the dynamic contribution from the hydro- 
mechanical portion of the servovalve. The load, consisting of the piston -in -cylinder 
actuator and the air valve, is represented by the expression in the third block. The 
mathematical expressions shown in figure 4 were developed in reference 8. The feedback 
transducer is represented by a constant and does not contribute any dynamics to the sys- 
tem over the frequency range of interest. 

A transfer function relating air valve position to input signal can be obtained from 
figure 4 by use of standard block diagram reduction techniques. The closed-loop transfer 
function for the system is 
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(All symbols are defined in appendix A. ) With the use of the constants for the system 
(listed in table I), equation (1) was evaluated over the frequency range of 1 to 200 hertz. 
The analytical result (fig. 5) indicates the existence of a resonance in the frequency 
range of interest. Later experimental tests verified the prediction of the model, and a 
resonance was encountered at approximately 150 hertz. This resonance was undesirable 
for two reasons: (1) valve damage could result because of a collision of the valve piston 
with its mechanical stops, and (2) for experimental testing, a constant valve stroke over 
the frequency range is desired. An improvement in system response was obtained through 
the use of a compensation network placed in the forward loop. The compensator char- 
acteristics in the frequency domain are represented by 


1 + 

G(jco) = m, 

i + J . q„ _ 

13 100 


( 2 ) 


Additional frequency response tests were run, and results showed that the compensation 
network did improve the system response. A comparison of valve response before and 
after compensation is shown in figure 6. The analytical results (dotted line) compare 
favorably with the actual system response. 


Pneumatic System 

The objective in the evaluation of the pneumatic system was to obtain large flow 
amplitudes from the distortion system over the frequency range of 1 to 200 hertz. A 
sketch of the pneumatic system appears in figure 7. The equations describing the dynam- 
ic behavior of the pneumatic system are listed in appendix B. These equations are 
written in terms of small perturbations about an operating point. The supply and dis- 
charge line dynamics are represented by transmission line equations developed in ref- 
erence 10. The transfer function relating jet nozzle pressure or flow (terminating 
orifice is choked) to valve area change is developed in appendix B and is represented by 
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The coefficients Cj, Cg, Cg, and C 4 , as defined in appendix B, are 



Cl = 1 + - K^a d 

(B14) 


C 2 = - K p3“d> 

(B15) 


C 3 = -R d (l + + K^ 3 a d ) 

(B16) 


C 4 = -R d R i( -l + Kp2 a i - K; 3 a d ) 

(B17) 

and 

Ri . 1 + ^ 
1 - “i K pl 

(B9) 


1 - a . 

R„- 

“d + l 

(B12) 

where 

erence 

a, a function of Mach number in the inlet or discharge line, 
10, is 

as developed in ref- 


o= 7M 
5 + M 2 

(B18) 


^ \ "O'/lS ( -(o.+t.)s\ 

jL)e d (l-Kje ' 1 / 

“d/ 


-(di+TjJs -(<r d+ T d )s 


+ C 3 e 


+ C 4 e 


"( CT i +T i +CT d +T d) s 


■(a.+TjJs 


The coefficients of the terms containing e represent the inlet line dynamics, 

- (°d +T d) s 

and the terms containing e represent the discharge line dynamics in equa- 

tion (B13). The constants and appearing in equations (B14) to (B17) repre- 
sent the valve sensitivity coefficients and are developed in appendix C. The term R. 
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is the inlet line reflection factor, and R d is the discharge line reflection factor. Con- 
sider the case where R i = 0, implying an impedance match for the inlet line, and the 
valve operating as a choked orifice (i. e. , = 1, = 0, K y = 1) . Under these con- 

ditions, the coefficients of the denominator of equation (B13) become 


C l = l 


+ a. 

l 


c 2 = ° 


(3) 


(4) 


C 3 " " R d^ + a i > 


C 


4 


= 0 


The system transfer function then reduces to 





2a d e 


~°d s 


*( a H +r H) S 

(1 + a d )(l + a.)|_l -R d e d d 


(5) 

(6) 


(7) 


The dynamic contributions are due to the discharge line alone and can be minimized by an 
impedance match in the discharge line. An impedance match is realized by allowing a? d 
to approach 1 (i.e. , R d —0). 

An evaluation of equation (B13) over the frequency range of 1 to 200 hertz for dif- 
ferent R^’s indicates how inlet line effects degrade system response. Figure 8 shows 
that large impedance mismatches in the inlet line (R^ — 1) result in attenuation at fre- 
quencies of 50 and 150 hertz. The degree of attenuation is controlled by the magnitude of 
the impedance mismatch: In later tests run on the system, results showed that the inlet 
line produces severe attenuations at frequencies of 43 and 129 hertz. The discrepancy 
between the analytical model and the actual system can be attributed in part to the Mach 
number chosen for use in the model. 

IMPEDANCE MATCHING 

The impedance of the reflecting surface must be made e^ual to the transmission line 
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impedance in order to provide an impedance match to eliminate the reflections. An 
orifice that has the proper impedance when placed in the supply line will eliminate re- 
flections. For an impedance match to exist in the supply line, R (eq. (B9)) must equal 
zero, which requires that 


oi. — - 

1 


K' 


( 8 ) 


Pi 


Substituting for K pl from equation (Cl 1), with y=1.4, results in 


a. = 


(2.4 - 2.8 K -1 + 2.8 


Pi 


K, 


tl, 


Pi 


1.4 K. 


(9) 


Pi 


where K^, defined by equation (C4), is rewritten using equation (Cl) and is rearranged 
to yield 


C d A = 


(0. 7143 - K pl ) 


ta 

3 sl/ 


.0.7143 ,; 7 2 


W T ta R 
P sl P ta g 


( 10 ) 


Equation (10) yields a value of effective -orifice area that will produce an impedance 

match in the supply line and is seen to be a function of the operating point. Figure 9 

shows the non -normalized response of jet nozzle pressure to valve area variation as a 

function of frequency for three different operating points. The zero-to-peak value of 

o 

valve area variation was 1. 6 square inches (10. 3 cm ) for all three cases. The open 
symbols denote the response without the impedance matching orifice. Supply line dy- 
namics produce antiresonances at approximately 43 and 129 hertz. With the use of 
equations (9) and (10) and a value of C^ = 0. 6 for the orifice, the area resulting in an 
impedance match for each operating point was determined. The required orifice area 
is included in figure 9. Later tests showed that the 4. 9 -square-inch- (31. 6 -cm -) area 
orifice produced the best dynamic results for all three cases. The solid symbols denote 
the improved response obtained with the orifice placed in the supply line and represent 
the response of the present system. The operating conditions shown in figure 9(a) result 
in a lower average flow but provide the largest flow variations because the valve is more 
effective in producing large flow variations for a given valve area change. The amplitude 
attenuation in the range of 10 to 60 hertz, shown by the solid symbols in figures 9(a) and 
(b), occurs because the orifice used is different from the computed value required for 
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an impedance match. The largest attenuation in the range of 10 to 60 hertz exists for 
the conditions of figure 9(a) because this case represents the largest mismatch. As the 
operating point changes (valve area increases), the degree of attenuation is reduced 
because the mismatch is reduced. In figure 9(c), an impedance match is achieved, and 
there is no attenuation in the range of 10 to 60 hertz . 

Other dynamic contributions to the responses shown in figure 9 come from the dis- 
charge line, the characteristics of which are shown in figure 10. The discharge line 
has resonances at 55 and 165 hertz, and the amplitude of these resonances is controlled 
by the magnitude of R d (eq. (B12)). The term R d is a function of the terminating 
orifice used. The discharge line resonance at 165 hertz appears in the response shown 
in figure 9. The broad resonance at 55 hertz does not show up as well because it is 
close to the supply line antiresonance at 43 hertz. 


CONCLUDING REMARKS 

The results presented herein define the problems encountered in developing an air 
distortion system. The analytical models for the servo and pneumatic system provided 
valuable assistance in defining the problem areas and aided in determining what modi- 
fications would help to improve the response. From experience with the present system 
and predictions of the analytical model, certain modifications will be incorporated into a 
second generation system. These modifications include close coupling of the valve and 
accumulator and reduction of losses through reshaping the valve ports and entrances of 
the pneumatic system. These changes are expected to result in a flatter response of jet 
pressure to valve area variation over the 0- to 200-hertz frequency range of interest. 

In addition, the pressure levels are expected to be higher because of the reduced losses. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 12, 1970, 

720-03. 


8 



APPENDIX A 


SYMBOLS 


A 

a 


. 2 2 
area, in. ; m 

speed of sound, ft/sec; 
m/sec 


C 1’ C 2 

c 3 ,c 4 


G(jw) 


discharge coefficient 

coefficients defined in 
eqs. (B14) to (B17) 

function of frequency 


g 

I 

K FB 


gravitational constant, 
ft/sec^; m/sec^ 

input signal, V 

gain of position trans- 
ducer, V/in. ; V/cm 



coefficients defined in 
eqs. (C3) and (C4) 
and (C9) to (Cll) 


K 


1 


variable gain of servo 
amplifier, A/V 


K 2 

K 3 


constant of servovalve, 
in. /(A)(sec); 
cm/(A)(sec) 

constant of servovalve, 
sec -1 


R, 


R. 


T 

W 


X 

a 


v 

p 

cr 

T 

CO 


k 4 

L 

M 

P 

R 


gain of load driven by 
servovalve, sec - * 

length, in . ; m 

Mach number 

p 

pressure, psi; N/m 

gas constant, ft/°R; 
j/(K)(kg) 


Subscripts: 

a,j 

d 

FB 

i 


coefficient defined in 
eq. (B12) 

coefficient defined in 
eq. (B9) 

Laplace transform 
operator 

absolute temperature, 
°R; K 

weight flow, lb/sec; 
kg/sec 

valve position, in. ; m 

coefficient defined in 
eq. (B18) 

ratio of specific heats 
(1.4 for air) 

damping ratio 

density, lb/ ft 3 ; kg/m 3 

transport time defined in 
eq. (B19), sec 

transport time defined in 
eq. (B20), sec 

frequency, sec - * 


station locations defined 
in fig. 7 

valve discharge side, 
downstream location 

position transducer 

valve inlet side 
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valve 


nl natural frequency of servovalve 

n2 natural frequency Of load 

o orifice 

p pressure sensitivity coefficient 
s static conditions 

t total or stagnation conditions 

u upstream location 


v 

1,2,3 station locations defined in fig. 7 
Superscripts: 

steady state 
perturbation 
time derivative 

associated with constants as de- 
fined 
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APPENDIX B 


PNEUMATIC SYSTEM DYNAMICS 

Assuming constant area and accumulator pressure, flow through the orifice is given 
by 





(Bl) 


where is defined by equation (Cll). 

The supply line dynamics are represented by equations developed in reference 10: 






+ 




The linearized equation representing the valve is 


W 


^v 


t2 


t3 


r =K vr + K » 2 f^ + K » 3 f 

w v Ay B t2 ^3 


(B2) 


(B3) 


(B4) 


where K^, and K p 3 are defined by equations (C9) to (Cll). Discharge line dy- 
namics are 






(B5) 

(B6) 


The terminating orifice that operates as a choked station is represented by 
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Hi=5i 

*j p ti 


(B7) 


Substituting equation (Bl) into (B2) and substituting the result into (B3) yields 


t2 

? t2 


= -a. 


-(a.+r.)s 
1 + R. e 11 


1 - Rje 


■(a i+ T i)s 


5 


(B8) 


where 


R, = 


1 “ “i 


K P1 


(B9) 


Substituting equation (B8) into (B4) and rearranging results in 


W. 


W 


v 


1 + K >i - R i^ - K p2 a i> e 
-(cr.+T.)s 


-(a.+r.)s 


1 - Rje 


- K ^ + K* ^ t3 
“ K VT“ +K p3f~ 
\ P t3 


(BIO) 


Substituting equation (B7) into (B6) to eliminate Wj/Wj and then substituting this result 
and (B7) into (B5) to eliminate P^/P^ and W./W. yields 


W. 


W, 


1 + V^d 5 

2 a. 


1-R d e 


-( CT d +T d) s 


P ti 


(Bll) 


where 


R d = 


1 - a c 
1 - a- 


(B12) 


Simultaneous solution of equations (B5), (B6) ? (BIO), and (Bll) to eliminate 
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/y — 

and W v /W v results in the final equation which represents the response of jet nozzle 
pressure to valve area variation. 


P tj 

r^i 

A^. 

\ 



1 - Rje 


-(CTj+T^S 


C 1 + C 2 e 


-( <r i +T i )s 


+ C 3 e 


■(a d+ r d )s 


+ C 4 e 


•(a i+ T i+ a d+ T d )s 


(B13) 


where 


C 1 = 1 + “i 1 ^ - “d^3 
C 2 = -1^(1 - - K^ d ) 

C 3 - ^d* 1 + K p2“i + ^3“d> 

C 4 = -W- 1 + ^>2“i - ^3“d> 

and a, a, and r, as defined in reference 10, are 

1-4M 
1 + 0.2M 2 


(B14) 

(B15) 

(B16) 

(B17) 

(B18) 


L 

or = — — (B19) 

1 + M 

T = — — (B20) 

1 - M 
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APPENDIX C 


VALVE SENSITIVITY COEFFICIENTS 


Compressible flow through an unchoked orifice is represented by 


W 



The following expression is obtained by linearization of equation (Cl): 


(Cl) 


where 


W _ A o , v P tu _ P sd 

P u f~ P d f” 

w A o p tu p sd 


V=(^ I + 


V. 


/p x (y-i)/r 

sd \ 

P tu/ 


1 - 


sd' 


,(y-i)/y" 




tu/ 


* sd \ (r - 1)/r 


J 


K pd =± < 

* y 


1 - 


(y - 1) \ tu; 


1 - 


sd' 


(y-l)/yl 


tu/ 


(C2) 


(C3) 


(C4) 


For consistency, equation (C2) should contain a total downstream pressure. The ex- 
pression for total pressure at a point is given by 


P td P sd + 


W 


2gp sd A d 


(C5) 
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The orifice loss assumed for this valve is taken care of by assuming that static pressure 
is constant over the mixing region downstream of the valve. Equation (C5) is valid for 
low Mach numbers in the pipe downstream of the valve. 

From the continuity equation, orifice flow is equal to downstream flow. From 
perfect gas law, 

P qd 

—2— = constant (C6) 


Substituting equation (C6) into (C5) to eliminate P sd and then linearizing the result 
yields 


P sd ^ r+1 ^ P sd“ P td P td ^ r+1 ^ P sd" P td W 


2y(P M - p c h) 


td *sd' W 


Substituting equation (C7) into (C2) and rearranging terms results in 


where 


w , P, P + . 

^■^r + K puf- + K pdf- 


(y+ i)f — ) - 1 


[< r+ l)-2^]^Mi-2yK^)l 
*td J 


} 

K pu “ V 


-i 


[(TH-D-ayK^j^-d^yKpa) 


(CIO) 
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(Cll) 


^pd 

|[(r + i)-2rK pt (]!j-(i-2 y K pd )j 

where K u and K pd are defined in equations (C'3) and (C4). For a choked orifice, 
equation Ijci) is 



W = 


and (C2) to (C4) become 


and (C9) to (Cll) become 



g y 



(y+i)/(y-l) 


C d A o P tu 



(C12) 


ry r>j 


_ o + r tu 

(C13) 

_A p tu 


J* 

¥ 

1! 

H - 1 

(C14) 

1! 

O 

(Cl 5) 

T— 1 

II 

(C16) 


(C17) 

o 

II 

.u 

M 

(Cl 8) 
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TABLE I. - CONSTANTS FOR SERVOSYSTEM 


Variable gain of servoamplifier, K^, A/V 
Servovalve constant, Kg, in./(A)(sec); em/(A)(sec) 

Servovalve constant, Kg, sec * 

Constant associated with load, K^, sec -1 

Position transducer, K FB , V/in. ; V/em 

Corner frequency of servoamplifier, co. , sec - '*' . . 
Natural frequency of servovalve, u; nl , sec . . , 

Natural frequency of load, sec -1 

Damping ratio of servovalve, ^ 

Damping ratio of load, 


623 ; 1580 
... 840 
7290 ; 7290 
80 ; 31.5 
. . 17 000 
. . . 4600 
. . . 5040 
... 0.4 
. . 0.186 



!. Flow line, supply to torus manifold 

2. Torus manifold 

3. Inlet duct 

4. Flow line, torus manifold to control valve 
.5. Control valve body 
6. Flow line, 


:ondary nozzle J 


rq 

lBiC-68-2313 


18 


(a) Installation of air distortion valves in test chamber. 
Figure 1. - Facility installation. 





(b) Array of distortion jets in inlet duct. 
Figure 1. - Continued. 


C-68-2315 







Amplitude ratio 
normalized 
to 1 hertz, 



Figure 3. - Air valve piston-cylinder assembly, 


Figure 4. - Dynamics of air valve position servo. 


2 4 6 8 10 20 40 60 80 100 

Frequency, Hz 

Figure 5, - Effect of gain on air valve response (analytical model). 


Amplitude ratio normalized 
to 1 hertz, 



Station 


a 
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Figure 7. - Pneumatic system for air valve. Inlet line: length, 66 inches 
(168 cm); outside diameter, 4 inches (10. 16 cm); discharge line (total of 
nine): length, 40 inches (101.6 cm); outside diameter, 1 inch (2. 54 cm). 
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Figure 8. - Effect of inlet line dynamics on system response for analytical model. 
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(c) For impedance matching, valve area of 6.3 square inches (40.65 crrrl 
and orifice area of 4. 9 square inches (31. 6 cmO. 

Figure 9. - Experimental results illustrating improved response obtained 
with impedance matching orifice in supply line. 



Figure 10. - Comparison of analytical and experimental results for discharge tubes. 
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